Superfluid density of the ±s-wave state for tlie iron-based su- 
perconductors 
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Abstract. - We study the superfluid density of the ±s-wave state of the minimal two band model for the Fe- 
based superconductors and its evolution with impurity concentration. We show that the impurity scattering 
of the strong coupling limit induces the selfenergy of a generic form ImLimp{(S)) ~ iY+ ipco beyond a criti- 
cal impurity concentration Timp > Tcrit- This form of Iml.imp{(£>) causes the temperature dependence of the 



superfluid density [Ps{T) — p.s(0)] 
impurity concentration r„„y, < Ten 



i —yT — (3r . Combining with the full gap behavior of Ps{T) for lower 
, the ±s-wave state produces a continuous evolution of AA,(r): exponen- 
tially flat with increasing impurity concentration that is consistent with the measurements of 
the Fe pnictide superconductors such as M-11 1 1 (M=La,Nd,Sm,Pr) and Ba-122 with various dopings, except 
LaFePO which shows AA,(r) o= r''^ at low temperatures by recent experiment. Our results also demonstrate 
that the density of states (DOS) measured by thermodynamic properties and the DOS measured by transport 
properties can in general be different. 
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Introduction. - The recent discovery of Fe-based super- 
conducting compounds [1-3], has greatly spuiTed the research 
of unconventional superconductors. With a discovery of new 
superconducting (SC) material, the most impelling question is 
to determine the SC gap symmetry. However, despite intensive 
experimental efforts, the pairing symmetry of this new class of 
SC materials is not completely settled. 

Tunneling spectroscopy of Ref. [4,5], photoemission spec- 
troscopy of Ref. [6], and NMR nuclear-spin-lattice relax- 
ation rate measurements [7-11] seem to indicate a nodal su- 
perconductors (SCs). On the other hand, tunneling spec- 
troscopy of Ref. [12], angle resolved photoemission spec- 
troscopy (ARPES) measurements of Ref. [13-15], specific heat 
measurement [16] all support a fully opened s-wave type gap. 
As to the penetration depth measurements, early measure- 
ments of M-1111 (M=Pr, Nd, Sm) [17-19], and (Ba,K)Fe2As2 
[20] showed an exponentially flat behavior at low tempera- 
tures supporting a s-wave type SC. However, recent reports of 
Ba(Fe,Co)2As2 [21,22] shows a power law of r2-2.5 ^j^^j jj^ore 
recently measurement of LaFePO [23] shows a near T-linear 
temperature dependence. 

On the theoretical side, it is almost agreed on that non- 
phononic fluctuation is the most probable pairing interaction 
and several theoretical proposals for the possible pairing sym- 



metries of Fe pnictide SCs were put forward [24-28]. However, 
recent advance of theoretical studies by several groups [29-32] 
have provided convincing evidences that the ±s-wave SC state 
is the most promising candidate for the true pairing state of Fe 
pnictide SCs. Therefore, we will consider only the ±s-wave 
SC state in this paper 

Having s-wave gaps on each bands but with opposite signs, 
the is-wave SC state would display standard features of a s- 
wave SC for various SC properties at its pure state, for ex- 
ample, the exponentially flat temperature dependence of the 
penetration depth, which is consistent with some of experi- 
ments [17-20]. Another hallmark of a s-wave pairing state is 
the Hebel-Slichter peak and exponential drop of the nuclear- 
spin-lattice relaxation rate l/T\, which is absolutely inconsis- 
tent with the experimental reports of all Fe pnictide SC com- 
pounds [7-11]. It was quickly pointed out by several groups 
[28-32] that the opposite signs of the order parameters (OP) in 
the is-wave SC state with a help of impurity scattering could 
reduce the Hebel-Slichter peak and also mitigate the exponen- 
tial temperature dependence into a power law temperature de- 
pendence. In particular, the Ref. [29] showed the special im- 
portance of the strong coupling (unitary limit) impurity scat- 
tering in the is-wave SC state. The source of strong coupling 
impurities is also naturally expected when there are detects or 
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vacancies on Fe sites, of which the concentration should be also 
very small (order of a few percents) unless the superconductiv- 
ity itself is substantially or completely destructed. 

To sketch briefly the previous work [29], impurity scattering 
of the strong coupling limit in the ±s-wave SC state induces 
an impurity bound state inside the SC gap, which is however 
off-centered from zero energy due to the generic property of 
the ±s-wave state, namely, the unequal sizes of the s-wave OPs 
with opposite signs on different bands. This impurity bound 
state drastically modifies the fully gapped DOS of the pure ±s- 
wave state into a V-shaped DOS as in a d-wave SC. This V- 
shaped DOS is not only consistent with the DOS directly mea- 
sured by tunneling spectroscopy [4, 5] but also conform with 
the isotropic gap measured by ARPES [ 1 3-1 5 1 ; a d-wave state 
would give a V-shaped DOS but at the same time should show 
a strong anisotropy of the gap in the ARPES experiments. Fur- 
thermore, it provides a clear physical explanation for the origin 
of the r"(a 3) power law of NMR 1 /Ti measured for all Fe 
pnictide SCs [7-11]. This is an excellent theoretical success 
in that a noble interplay between the strong coupling impuri- 
ties and the ±s-wave gaps provides coherent resolutions to the 
several conflicting experimental observations with only one ad- 
justable parameter, i.e., the impurity concentration. 

In this paper, we extend the previous work [29] to study the 
superfluid density and penetration depth in order to complete 
the consistency test of the ±s-wave state as the ground state 
of Fe pnictide SCs, except LaFePO [23]. The main results 
of this paper is that the temperature dependence of X{T) of 
the ±s-wave state evolves systematically with increasing impu- 
rity concentration: from a exponentially flat behavior for pure 
sample to a behavior at a critical impurity concentration, 
and finally to a one beyond the critical impurity concen- 
tration. These temperature dependencies are in agreement with 
file data of M-11 1 1 (M=Pr, Nd, Sm) [17-19] and (Ba,K)Fe2As2 
[20] (exponentially flat), and Ba(Fe,Co)2As2 [21,22] (power 
law oz T^^^-^). However, recent report of AX{T) oc ji-^ for 
LaFePO [23] is not compatible with our scenario. Possible ex- 
planation for this exception is that this compound may have a 
different pairing symmetry than the other Fe pnictide SC com- 
pounds [28]. We will not pursue this issue in the current paper. 

Model and Formalism. - For the minimal two band 

model of the its-wave pairing state, we assume two s-wave 
OPs Aft and Ag on the two representative bands of the Fe pnic- 
tide materials: one hole band around F point and one elec- 
tron band around M point in the reduced BriUouin Zone. A/, 
and Ag have opposite signs and their magnitude are different in 
general. Impurity scattering will renormalize the energy (co) 
and the OPs (A^^e) through the selfenergy corrections: nor- 
mal selfenergy and anomalous selfenergy, respectively. These 
impurity induced selfenergies are calculated by the T -matrix 
method [33, 34], suitably generalized for the ±s-wave pairing 
model [29], as foUows. 
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ni. 



imp 
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Ah,e = Ah,e + '^l{(a„)+I.l{(a„), 



(1) 
(2) 



where a)„ = TK{2n + 1) is the Matsubara frequency, ni,„p the 
impurity concentration, and Ntot = Nh{Q) +Ne{Q) is the total 
DOS. The T -matrices t"'^ are the PauU matrices x"'^ compo- 
nents in the Nambu space and defined as foUows. 
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where c = cotSo is a convenient measure of scattering strength, 
with c = for the unitary limit and c > 1 for the Born limit 
scattering. (...) denotes the Fermi surface average. 

With the self-consistently calculated (h„ (Eq.(l)) and A/, g 
(Eq.(2)), we can calculate all physical quantities of the model. 
For small amount of impurity concentration (F/Ag <C 1) of the 
strong scatterers (c = 0), which we believe the case of the rel- 
evant experiments, the Tc suppression is small [29] and the 
renormalization of the OPs A^ ^ is also marginal. However, the 
renormalization of &>„ is strong even with a small amount of 
unitary scatterers and develops an off-centered impurity bound 
state inside the SC gap. This unusual low energy impurity band 
modifies the fully gapped DOS of the pure state into a V-shape 
DOS and becomes the origin of the power law (°= T^) of the 
NMR relaxation rate l/7i of Fe pnitide SCs [7-11]. In this pa- 
per, we study how this modified DOS of the its-wave state with 
impurities affect the superfluid density and penetration depth. 

Avoiding unnecessary repeat, we begin with showing pre- 
vious calculations of the DOS Ntot{^) and the correspond- 
ing normal selfenergy Imll^gf (co) = /mE^ -t /mE^ for various 
impurity concentrations in Fig. 1(a) and Fig. 1(b). These re- 
sults were calculated with the parameters |Aej/|A/,| = 2.5 with 
Nh (0) /Ng (0) = 2.6, which were obtained with the realistic band 
structure parameters [29]. The physical results of the present 
paper, however, are independent of the specific choice of the- 
ses parameters. All energy scales are normalized by \Ag\ in 
this paper. Fig. 1(a) shows the total DOS with different impu- 
rity concentrations T/Ae = 0.0,0.01,0.04,0.08 of the unitary 
scatterer (c=0), and Fig. 1(b) shows the corresponding impurity 
induced selfenergy /mZ^^j ((o) = Iml!^ + ImLg. 

Figure 1(a) shows systematic evolution of how the fully 
opened gap of pure state is filled with impurity states; the pat- 
tern of filUng is very unusual and the F/Ag = 0.04 case displays 
a perfect V-shape DOS down to zero energy similar to a pure 
d-wave SC DOS. The origin of this behavior is easily seen in 
Fig. 1(b); the impurity bound state is never formed at zero en- 
ergy but away from it because of the incomplete cancellation 
of [GI + GI] in D (Eq.(5)), the denominator of the T -matrices 
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T°^'. Were the sizes of A'/,|A/,| and Ne\Ai.\ equal, the FS av- 
eraged anomalous Green's functions [Gj^ + Gl] would vanish 
and then the impurity bound state would form at zero energy; 
this condition is intrinsically satisfied in the case of d-wave SC 
state. In generic ±s-wave state, therefore, the full gap around 
i£> = is protected until this off-centered impurity band spills 
over to the zero energy with increasing impurity concentration. 
When it touches the zero energy limit, the ±s-gap state will 
thermodynamically behave like a pure d-wave SC, and this hap- 
pens with the critical impurity concentration Ycrit {— 0.04Ae for 
our specific model parameters). Increasing impurity concen- 
tration beyond Ycrit, the DOS still keeps the V-shape but now 
Ntot{^ = 0) obtains a finite value (see the Y — O.OSAg case in 
Fig. 1(a)). 

Superfluid density and penetration depth. - With the 
above calculated (b„ and A/,.e, we can study the temperature 
dependence of the superfluid density along the Fe-plane, with 
the following kernel. 



^(7^) = - L Na2%TY^ 



a=h,e 



„Re 



(a2+A2)V2 



(8) 



The above kernel is directly proportional to the superfluid 
density Ps{T) and l/Xj^{T) in the London limit. We use this for- 
mula for numerical calculations. Before we show the numerical 
results, we would like to gain analytic understanding about the 
temperature dependence of the superfluid density p.v(r) of the 
±s-wave state when the DOS is modified by impurity scatter- 
ing as shown in Fig. 1 . To focus on the temperature dependence, 
we separate only the temperature dependent part from Eq.(8) as 



5K{T) = -— ^ Njvl. r^co/(co)Re- 



■(a)2-A2)3/2 



(9) 

where /(co) is Fermi-Dirac distribution function and cb and 
Aa(a)) are the real frequency quantities obtained from co,, and 
Aa(a)„) by analytic continuation (co,, — ^ co + /r|). As explained 
before, for small impurity concentrations, A„ w A„ and A„ 
doesn't develop any peculiar low energy structure, therefore we 
focus on the behavior of CO = co + E^gj(co) at low energies. As 
can be seen in Fig. 1 (b), for r,-,„p > Ycrit (^crit = 0.04Ap), we can 
approximate CO « aco + /pco + iy... for low frequencies. Here 
Y is the typical impurity induced constant damping. Unusual 
part is the term /pco which gives a linear-in-co damping that is 
a unique feature of the ±s-gap state with impurity scattering 
of a strong coupling limit (unitary limit). With this form of CO, 
we can extract the temperature dependence of dK{T) using ex- 
tended Sommerfeld expansion [35] at low temperatures. After 
straightforward calculation, we obtain the following result. 
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Fig. 1: (Color online) (a) Normalized DOS Ntoi{_&) for different im- 
purity concentrations, T/\ = 0.0,0.01,0.04,0.08. Thin dotted line is 
of the pure state for comparison and other lines are offset for clarity 
(the zero baselines of the offset are marked by the narrow horizontal 
bars of the corresponding colors), (b) Impurity induced self energies 
/raZ^^,(co) = /raZ^ + /mZ^ with the same parameters as in (a). These 
curves are not offset. 



As expected from a power counting, the constant damping term 
y produces decrease and the linear-in-co damping term pco 
produces decrease of superfluid density f>s{T). The de- 
crease of pi(r) due to a constant damping is the well known 
coiTection due to the strong coupling impurity scattering in the 
d-wave high-r^- cuprates [36]. What is unusual is the cor- 
rection due to the linear-in-co damping. This linear-in-co damp- 
ing is the unique feature of the ±s-wave SC with strong cou- 
pling (unitary limit) impurity scattering and it modifies other- 
wise fully opened DOS into a V-shape DOS. This V-shape DOS 
results in the same thermodynamic behaviors as in a d-wave 
SC such as in specific heat and NMR \IT\ relaxation rate, etc. 
However, when transport property like the superfluidity density 
is calculated, specific vertex enters combined with the coherent 
factors and a naive picture based on the DOS fails. 

With this result we have now a complete evolution of the 
temperature dependence of the superfluid density (T) for the 
lbs-wave SC with strong coupling impurities. For a clean sam- 
ple {Yiinp < Ycrit), it will show an exponentially flat dependence 
at low temperatures. With increasing impurity concentration. 
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Fig. 2: (Color online) (a) Superfluid density Ps{T) for different 
impurity concentrations, T/Ag = 0.0,0.01,0.04,0.08, normalized by 
Pj(r = 0) of the pure state, (b) Corresponding penetration depth 
A^(r) = X{T) — X{T = 0) normalized by Xq{T = 0) of the pure state. 
The data for F/A^ = 0.04,0.08 are offset for clarity. The power law 
lines (sold black lines) of and are shown for comparison. 



this flat temperature region shrinks. At the critical impurity 
concentration (Ti^p = Tcrir), it will display the behavior 
And finally for r„„p > r„it, the behavior will smoothly 
dominate over the dependence. Therefore, for a finite win- 
dow of impurity concentrations of F,-,,,;, > Ta-it, experiments 
would display a power law oc with 2 < a < 3. This evolu- 
tion of temperature dependence with impurity concentration is 
consistent with flie data of M- 11 11 (M=Pr, Nd, Sm) [17-19] 
(flat), (Ba,K)Fe2As2 [20] (flat), and Ba(Fe,Co)2As2 [21,22] 



(oc r^-^-^), except for the data of LaFePO [23] (°<: T' ^. 

In Fig. 2, we show the numerical results of K{T) of 
Eq.(8) using the same parameters as in Fig.l; K{T) is the 
same as the superfluid density Ps{T) and proportional to 
1 /Xj^{T). For the temperature dependence of the gaps A/, 
we used a phenomenological formula, Af,^e{T) = Af,,e{T — 
0)tanh(1.74-y/rc/r — 1). The values of Ah,e/Tc determine the 
high temperature behavior below T^. We chose lAjjTc = 3.0 
and |Ag|/|A/,| = 2.5 for our calculations. 

Figure 2(a) shows the superfluid density p.v(r) over 



Q < T < Tc for different impurity concentrations, F/Ag = 
0.0,0.01,0.04,0.08, normalized by p°(r =0) of the pure state. 
We warn the fact that the temperature dependence of Pi (r) at 
high temperatures (say, 0.5 < T /Tc < 1) does not reflect the 
low energy DOS nor the impurity effects but it is mainly de- 
termined by the temperature dependence of A/, ^(r). Ps{T) of 
the pure case (F/Ag = 0.0) shows the exponentially flat tem- 
perature dependence at low temperatures, consistent with the 
fully gapped DOS (see Fig. 1(a)) of the ±s-wave state. With 
increased impurity concentration, F/Ag = 0.01, this full gap 
behavior still persists. When impurity concentration becomes 
critical, i.e., T/Ag ~ 0.04, the DOS become a perfect V-shape 
down to zero frequency and the exponentially flat region of 
Ps{T) disappears. With further increase of impurity concen- 
tration, for r/Ae ~ 0.08, the DOS is not only a V-shape but 
it also obtains a finite DOS at zero frequency N{Q) = const., 
and p.s(r) becomes even more smooth. To see more clearly 
the power laws of the low temperature behavior, we plot, in 
Fig. 2(b), the temperature dependent part of the penetration 
depth AA.(r) = X{T) - X{T = 0) normalized by Xo{T = 0) of 
the pure state. The power law fittings (black solid lines in 
Fig. 2(b)) clearly confirm that AA,(r) °^ at the critical im- 
purity concentration F;,„p = Tdt = 0.04Af, and AX{T) °^ for 
higher concentration F,„,p = 0.08Af,, as we have shown with an 
analytic analysis of Eq.(9). 

Thermodynamic DOS and transport DOS. - Switching 
gears, we would like to draw attention to the more general 
aspect of DOS, which was clearly demonstrated in this pa- 
per; namely the fact that the DOS measured by thermodynamic 
properties and the DOS measured by transport properties can 
in general be different. The DOS of the its-wave state with 
a critical impurity concentration looks a V-shape as shown in 
Fig. 1(a). Here the DOS is the distribution of the quasiparti- 
cle energy eigenstates and this distribution can be measured 
by thermodynamic properties such as specific heat, photoemis- 
sion, and spin-lattice -relaxation rate l/Ti, etc. Therefore, this 
V-shape DOS is thermodynamically indistinguishable from the 
V-shape DOS of a pure d-wave SC state. As mentioned, how- 
ever, angle-resolved-photoemission spectroscopy can still dis- 
tinguish the presence or absence of the angular (or Fermi sur- 
face) anisotropy of the quasiparticle DOS [13-15]. 

On the other hand, the superfluidity density is a transport 
quantity defined as 



J{r) = ~-K{T)A{r) 

= -—p{T)A{r) 
mc 



A{r) 



(11) 



in the local limit (London limit). In the non-interacting case (no 
further quasiparticle scattering after diagonalization of pairing 
interactions), the BCS theory give a simple formula of the ker- 
nel to the vector potential A as [37] 



1-2 d03 
A 



CO 



3(0/ ''(co2-A2)1/2 



(12) 



p-4 



Superfluid density of the is-wave state for the iron-based superconductors 



The second term is apparently proportional to the thermo- 
dynamic DOS Re f^^2 a2)'/2 when 6) = + the gen- 
eral expression of the kernel (Eq.(9)) can be converted into 
Eq.(12). However, once the quasiparticle obtains a finite life- 
time due to interactions, Eq.(9) and Eq.(12) can not be inter- 
changeable and we have to use Eq.(9) to calculate the super- 
fluid density and penetration depth. Even for the simplest case 
of a constant damping O) = (H + t^, using Eq.(12) would give 
5X{T) oc —const., but the correct temperature dependence due 
to a constant damping is Sk{T) °^ —jT^ as well known in the 
case of high-Tc cuprates 1 36] and reproduced in the current pa- 
per When CO — co + ZpCO, the difference is even more interesting; 
Eq.(12) predicts dX{T) oc — pr but the correct behavior from 
Eq.(9) produces dX{T) °^ —^T^. We emphasize that this dif- 
ference of the DOS depending on the different measurements 
should be carefully considered for interpretations of various ex- 
periments of Fe pnitide SCs. 

Conclusion. - We have calculated the superfluid density 
Ps{T) and penetration depth X{T) of the ±s-wave state with 
impurities of the strong coupling hmit. We showed, both by 
analytic analysis and by numerical calculations, that the tem- 
perature dependence of ps{T) and X,(r) at low temperatures 
continuously evolves in a sequence of the forms: exponentially 
flat — — with increase of impurity concentration. 
This result consistently explain the experimental data of M- 
1111 (M=Pr, Nd, Sm) [17-19] (flat), (Ba,K)Fe2As2 [20] (flat), 
and Ba(Fe,Co)2As2 [21,22] (oc r^-^ S) xhe near T-linear be- 
havior of LaFePO [23] data can not be explained by our theory. 
If the data of LaFePO [23] is indeed confirmed, it indicates that 
LaFePO may have a different pairing symmetry than the other 
Fe pnictide SCs. Finally, we noticed that Vorontsov et al. [38], 
have recently studied the same problem as in this paper, and 
obtained a similar result PsiT) «; for high concentration of 
impurities but obtained a different result Ps{T) o= T^-^ for the 
critical impurity concentration. We think that the origin of this 
difference mainly arises from the fact that Vorontsov et al. [38] 
studied a weak coupling impurity scattering and we studied the 
strong coupUng impurity scattering. As a consequence, the def- 
inition of the critical impurity concentration and the manner 
that the SC state becomes gapless with impurities is different 
in each studies. 
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